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ABSTRACT 


The propagation of grazing incidence plane waves along a finite impedance boundary is 
investigated. A solution of the semi-infinite problem, where a harmonic motion, parallel to the 
boundary, is imposed along a line perpendicular to the boundary, is obtained. This solution consists 
of quasiplane waves, waves moving parallel to the boundary with amplitude and phase variations 
perpendicular to the boundary. Several approximations to the full solution are considered. 


INTRODUCTION 


Mathematical modeling of the propagation and reflection of harmonic plane waves above a 
finite impedance plane surface is a fundamental topic in acoustics. In the case where the angle 
between the normal to the wavefront and the surface is not zero an analytic solution is very easy to 
obtain. This solution consists of the incident plane wave propagating toward the surface plus a 
reflected plane wave propagating away from the surface at the same magnitude of the angle between its 
normal and the surface as the incident wave. The amplitude of the reflected wave is given by a 
reflection coefficient that is expressed in terms of incident angle and the specific impedance of the 
surface. However at zero incident angle (the wave normal parallel to the surface), complete cancellation 
of the incident and reflect waves occurs in this model and a zero solution results. Most acoustic texts 
claim that this situation is not possible [1-3]. 

McAninch[4] recently has investigated a related situation where a plane wave source is 
generating waves that would move parallel to a surface if its impedance was infinite but where the 
surface impedance isnot infinite quasiplane waves result McAninch's investigation, however, uses the 
parabolic approximation where only waves traveling in one direction are allowed. This paper 
approaches the same problem without the assumption of parabolic approximation. 


FORMULATION OF THE PROBLEM 


The governing acoustic wave equation for harmonic waves can be put in the form 

( V 2 + k 2 ) <J>=0 (1) 

where the time dependent part of the potential, e- i(0t , has been separated from the spatial part of the 
potential, 4>(x,y). When an impedance boundary exists, the solution of equation (1) must also satisfy 
the boundary condition 

4> y + Y<1> = 0 (2) 

on y=0. Here the subscript y indicates a partial derivative with respect to y. 
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For a unique solution, some extra constraints must be introduced. One is to assume that <|) will 
not be affected by the ground impedance as y approaches infinity, the second is that the acoustic 
pressure at x=0 is given by 


<K0,y) = 1 


(3) 


EXACT SOLUTION 


We assume that the solution of (1) has the form of 

([) = e ikx + f(x,y) (4) 

where e ikx can be considered as a solution without the boundary condition given by (2). Substituting 
(4) into (1),(2) and (3) we get a new governing equation and set of boundary conditions 

( V 2 + k 2 ) f = 0 (5) 

f y (x,0) + yf (x,0) = - ye 1 k x (6) 

f (0,y) = 0 (7) 

and 

lim f (x,y) = 0 (8) 

y— 

Equation (8) results from the first uniqueness condition listed above. 

The sine transform, 

oo 

F (X,y) = J f (x,y) sin (kx) dx (9) 

o 

is equivalent to the Fourier transform of an even function and will be applied here. The inverse 
transform is given by 


f (x,y) = ~ J F (x,y) sin (kx) dk 
0 


Applying (9) to (5) and (6) we have 


and 


3 2 F 

3y 2 


( k 2 - k 2 ) F = 0 


F (k, 0) + yF (k,0) = - y— r^— r 

7 1 -k 2 


( 10 ) 


(ID 

( 12 ) 
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Here, it is assumed that k is a complex number with a very small positive imaginary part. 
Solving (11) and making use of the given boundary conditions, yields 


F (X,y) = A (X,0) e my 

(13) 

A(K 0 )~- 

( X. - k ) ( y- m ) 

(14) 

where m = V(X 2 -k 2 ) . Since the solution is required to remain finite, ReV(X 2 -k 2 ) >0. Substituting (13) 
and (14) into (10) yields the inverse transform of F(X,y) as 

oo 

f (x,y) = - -rp f — — — e my sin (X,x) dX. 

" i(X 2 - k 2 )(y-m) 

(15) 

For convenience, substitute the identity 


. J Xx i Xx 

sin (Xx) = — — ^ 

2 1 

into (15), yielding 

(16) 

f( x,y ) = - 1 X (I ,- I 2 ) 

(17) 

where 


oo 

r 1 ’ my 

•l = J“T! eiX%dX 

o (*> ' k )(7-m) 

(18) 

oo 

r 'l _ m y 

>2 = [—5-“ 

o (^ - k )(y-m) 

(19) 

In order to evaluate the above two integrals, introduce the complex variable A = 

X+is and define the 


contour integrals 



Ae 


-My 


( A 2 - k 2 ) ( y- M) 


Ae 


-My 


( A 2 - k 2 ) ( y- M) 


e iAx dA 
e 1 Ax dA 


( 20 ) 


( 21 ) 


First evaluate the integral Iq where the contour is shown in Figure 1 along with the branch lines 
which extend from the imaginary axis to the points A = ± k. The value of this contour integral is 
determined by the residue within the contour. Writing 


Y= a + i |3 


( 22 ) 
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it is clear a pole exists within the contour only when a > 0 (since Re^A i 2 - k 2 > 0, y - VA 2 -k 2 =0 only 
when Re(y) = a > 0), and this pole is at A = \(k 2 + y 2 ). It is easy to determine the residue at this pole 


to be 


Res 


(Vk 2 +f )=-£ 


yy +i 7 k 2 +-y ? x 


and 


*1 " *C1 ' ^C2 ' ^C3 ' *C4 ' ^C5 ' *C6 + 


2 7t i Res ( -J k 2 + y 2 ) a > 0 
0 a < 0 


when R — > co, l C2 will vanish, while 

OO 

T C3 = ' J 


-i 7 s 2 + k 2 


y e sx ds 


o ( s 2 + k 2 ) ( y- i 7 s 2 + k 2 


(23) 


(24) 


(25) 


C4 


K. 

"J 


-i 7k 2 - X 2 y j Xx 


2 ' 2 ^-- iTi 2 ^ 2 


e' dX 


o ( X - k 6 ) ( y- i V k' - X ) 


t - . ILL J k x 
X5 y e 


I 


C6 


K 

-J' 


i 


‘ 2 ' 2 ivV 


e 


o (X. -k ) (y+iVk -X. ) 

Substituting the above integrals into (24), yields 


K 




.2 \ 2 
k - A, y 


-i 7k 


k 2 -\ 2 y 

6 ^ e‘ Xx dX 


J(?i 2 -k 2 ) [(y+i7k 2 -^ 2 ) ( y- i Vk 2 - X. 2 ) 


i 7 s 2 + k 2 


oo 

I — 

o ( s 2 + k 2 ) ( y- i v s 2 + k 2 


e' s x ds + 2-ZL e ikx 


+ 


i 2 n - yy J -/kV^ 

Y 

0 


a > 0 
a < 0 


(26) 


(27) 


(28) 


(29) 
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The value of integral I 2 is much easier to evaluate. We choose a contour in the fourth quadrant, since 
there is no pole within the contour. I 2 can then be written as 



s e 


ij 


s 2 +k 2 y 


( s 2 +k 2 ) ( Y+iTs 2 


+ ) 


e' sx ds 


(30) 


Here it must be recalled that there is a branch line along the imaginary axis. Subtracting I, from I, 
yields 



- i V s 2 + 


k 2 y 


i 777 


k 2 y 


00 

f s 

o ( s 2 + k 2 ) [ ( y- i V s 2 + k 2 ) ( y+ i 7 s 2 + k 2 ) 


>e sx ds 


+ 




i 271 -yy i •/ k 2 +y 2 x 

y c c 
0 


a > 0 
a < 0 


(31) 


By substitution of A.=Vk 2 -t 2 and s=-W k 2 -t 2 , the above two integrals can be combined into one. The 
final result is 


where 


and 


f(x,y) 


- e‘ k x + P - K 
-e ikx -K 


a > 0 
a < 0 


P = 2 e 


- yy 



X 


(32) 


(33) 


K= 2jr f ( t Cos (ty) - YSin (ty) ) ^ V k 2 . t z x dt 
71 g t(f+t 2 ) 


(34) 


P is called the surface wave, and it both decays with increasing height y, and also decays with the 

distance x due to the imaginary part of Y 
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ASYMPTOTIC EXPANSION VALID FOR SMALL RECEIVER HEIGHTS 


1. Soft boundary case 

Integral K can be asymptotically evaluated for large x using the saddle-point method. This 
method is discussed by Morse and Feshbach[5] and will not be discussed here. Actually we can use 
some of conclusions from Wenzel [6] since we have the same factor Vk 2 -t 2 as occurred there. 

The steepest-descent path has been shown in [6] to be given by 

T = t + i s (35) 

where 

(36) 


and t > 0. Again using the residue theorem, integral K can be transformed into the integral L. Note 
that i y is not in the region of concern since |3>0. Thus 



K = <! 


f P + L 


- i ye D 

- i Y £ D 


(37) 


where D is the region between positive real axis and curve s=-t (1+t 2 / k 2 ) 1/2 , P is the surface wave 
given in (33) and 


2 Y f T Cos (Ty) - y Sin (Ty) x dT 

71 sit (f+T 2 )T 


(38) 


Substituting (37) into (32), we have 



P + L 
L 


- i y e T 

- i Y <£ r 


(39) 


The region T in y plane is bounded by the curve p > 0, 0 < a=p(l+ (3 2 / k 2 ) 1/2 . The region T is 
called the surface wave region in the far field (shown in Figure 2), which is same as that of reference 
[6]. When ye T, we can easily show ReV k 2 + y 2 > k,that means, if the surface wave exists,its 
propagation speed is less than the speed of sound in free space. It is also found that Im "v k 2 + y 2 
has a close relationship to the quantity (a P / k) so that a large imaginary part of y and low frequency 
of the source can make the surface wave decay very quickly. 

Expanding each expression in (38) around the saddle point T=0 and integrating each term, we 
get 
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( 40 ) 


=i /jO 

Y V in 


i lc T 3 y 3 9 

Jkx<!(l-Yy) + -LL (1- Yy+ 1JL.JL^) + 0(X- 2 ) 
I x 


This asymptotic expansion is not uniformly valid. The conditions for its validity are 


and 



(41) 

(42) 


Extremely small l Y I will not satisfy conditions (41) and (42), so another asymptotic method has to be 
developed. 

The total solution under the condition of large x can be obtained by substituting (39) into (4) 
yielding 


f P-L 



y € r 

Y e r 


(43) 


where P and L are given in (33) and (40). If we neglect the surface wave, we can get an explicit 
equation for the wave above ground in the far field as 


(j> = 


lYl 


7 ? 


((l-ay) 2 + (Py) 2 ) 


i(kx+^- +0) 
_ 4 


(44) 


where 0 = - arctan (p y / (1 - a y)). Furthermore, if the receiver is on the ground, the above 
expression can be written as 


20 Log <J>= 20 Log a - 10 Log x (45) 

with a = (l/lY0(2k/7i) 1/2 . This result shows that the acoustic pressure level drops lOdB when the 
distance increases 10 times or 3 dB per doubling of distance. 

2. Hard boundary case 

As mentioned before, the asymptotic expansion given in (40) is not uniformly valid in y, with 
the method failing for small lyl. An alternative method is developed in this section which is valid in 
the small lyl case. The method is almost the same as that used in evaluating L except the factor 1 / 
(Y 2 + T 2 ) in (38) will not be expanded. After changing variables (38) becomes 


L = 


Y ( 1 - YY ) 

7t i 


Making use of the formula 



A kx 


I< 


X “jp 

2 i k 



(46) 
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( 47 ) 


yry e~ z2 erfc ( -iz ) Im (z) > 0 

-~-t z ( erfc ( -iz ) - 2 ) lm(z)<0 

2 1 z 

and neglecting the terms of order y 2 , yields 

-e ikx [l+y(-l r /IT _ y)] Im(-?-)>0 

Jti v 2 k yi 

(48) 

. e i kx [i + y(^_ /IT -y)-2(l-yy)] Im(-|-)<0 

V TC V 2 K yj 


The conditions Im (y/ Vi ) > 0 and Im (y/ Vi ) < 0 can be identified as a < (3 and a > (3 respectively, 
a = (3 is the line which divides these two regions in y plane. This is exactly the bounding curve a = [3 
(1 + P 2 /k 2 ) _1/2 obtained previously provided that lyl — » 0. Recognizing this relation, we substitute (48) 
into (43) and rewriting surface wave approximately as P = 2 (1 - y y) e ik \ finally get the total field 
expression as 





the condition for the validity of the above expansion is 

lY|y«l 

and 

lir | V ! « 1 

although x can't be small because of the nature of the saddle point method. 


(49) 


(50) 


(51) 


ANOTHER ASYMPTOTIC EXPANSION VALID FOR LARGE RECEIVER HEIGHT 


The asymptotic expansions obtained above have their limitations in application. For example, 
they require the receiver's location to be near the ground. In this section we will derive a asymptotic 
expansion which is valid for large R = Vx 2 +y 2 (except for small y) . The idea is similar to that of 
Chien and Soroka [7]. 


By using the identity 

f(x,y) 


sin (Xx)=(e‘^ x -e' i ^ x )/2i and transformation 
_ y f Tan ( z ) e ik ^ y Cos(z)+x 

in J y + i k Cos ( z ) 


X=k sin (z), (15) becomes 
s,n<l)) dz (52) 
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The contour C is shown in Figure 3. In order to get an expansion in terms of the variable R = Vx 2 + y 2 , 
we transform the Cartesian coordinate system into the polar coordinate system by 

x = R Sin 0 (53) 

and 

y = R Cos 0 (54) 

Substituting into (52), yields 


f (x,y) = -X- f — ^ 

1 71 J Y+ 1 


Tan (z) 


y + i Cos ( z ) 


ikRCos(z-G) j 

e dz 


(55) 


the saddle point for the function i k R Cos(z-G) is at z = 0 and the path of steepest descent is found to 
be given by 


Cos ( u - 0 ) Cosh ( v ) = 1 (56) 

where z = u + i v, by considering Im (i k R Cos (z - 0))=lm (i k R). This path, denoted as C' is 
shown in Figure (3). Deforming contour C into C', adding the possible poles (Cos z = i y / k ), we 
have 


f (x,y) = Q + H ( - Re ( 1 - ^ Cbs 0 - 

Jv 



Sin 0 ) ) P 


(57) 


where Q is defined by (54) but with the contour C changed to C', H is Heavyside step function and P 
is the surface wave given in (33). The condition for the existence of pole is explained in reference [7], 
and will not be repeated here. In the limit of 0 approaching 7t/2 the condition for the existence of the 
pole in the present case is equivalent to the condition for the existence of the pole in (39). 

Q can be evaluated asymptotically with a method similar to that used in evaluating L, i.e. to 
expand each term around the saddle point 0 and then integrate them with suitable transformation of the 
variable . The result is 


q= /24— — — 

V 1 7t k R i ( y+ i 


yTan 0 ikR 
c 


1 +T 


1 


ikR 


1 + 


Cos 0 ) 
i k Cos 0 


i k 


2 ( y+ i k Cos 0 ) ( y + i k Cos 0 ) Cos 0 


The conditions for the validity of the above expansion are 

kR » 1 


and 


Cos 2 0 , ^ 


y + i Cos 0 R Cos 0 


« 1 


(58) 


(59) 

(60) 
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It is clear that 0 cannot be too close to Till. This limitation is complimentary to the asymptotic 
expansions obtained previously (for small y). Substituting into (4) yields 


iJ)=e ikx +Q + H(-Re(l - Cos 0 - 1 + “T Sin 0 )) P (61) 

In the limit R->oo Q and P will vanish , with the result that only the plane wave term remains. 

Equation (38) can be evaluated accurately by numerical methods as well as by asymptotic 
expansions. Calculations show that the results match quite well when y is small. Figure 4 gives the 
amplitude of acoustic pressure on the ground versus the distance to the receiver obtained by numerical 
integration and from (44). Figures 5a, b and c show the amplitude of acoustic pressure versus the 
receiver height for several receiver locations as obtained from the asymptotic expansions, (40) and 
(61). These figures are similar to the results obtained by McAninch [4], 


CONCLUSIONS 


The acoustic field of a plane wave at grazing incident to a finite impedance has been 
theoretically investigated. Exact numerical and asymptotic expansions are developed, which are very 
similar to those found by Wenzel [6] for a point source and by McAninch f4j using the parabolic 
approximation to the wave equation. When y is small, the incident wave is indeed canceled, but the 
result is not zero due to the existence of a surface wave and the wave denoted as L. Near the ground, 
the acoustic pressure decays as x' 1/2 (assuming the surface wave is neglected). The asymptotic 
expansion for large distance R shows that the acoustic pressure decays as R 1/2 when R— and 
when the receiver is not close to the surface only incident wave exists. 
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Figure 3. The u-v plane, showing the integration contour C and the steepest descent path C'. 
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Figure 4. Comparison between a numerical result and the asymptotic result, (44), for the parameters, 
k = 18.313 and y = 6.935 + i 19.015. ’ ‘ 
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